'Hole-digging' in ensembles of tunneling Molecular Magnets 
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The nuclear spin-mediated quantum relaxation of ensembles of tunneling magnetic molecules 
causes a 'hole' to appear in the distribution of internal fields in the system. The form of this hole, 
and its time evolution, are studied using Monte Carlo simulations. It is shown that the line-shape 
of the tunneling hole in a weakly polarised sample must have a Lorentzian lineshape- the short-time 
half- width £ D in all experiments done so far should be ~ Eo, the half- width of the nuclear spin 
multiplet. After a time t , the single molecule tunneling relaxation time, the hole width begins to 
increase rapidly. In initially polarised samples the disintegration of resonant tunneling surfaces is 
found to be very fast. 

PACS numbers: 



In the last few years many experiments have been 
done on the tunneling relaxation of magnetic molecules 
jH [J [1 @ These couple to each other, via dipole 
and weak exchange interactions, and also to phonon 
and nuclear [j| environments. The most complete results 
have appeared on the F e-8 system, where measurements 
of the relaxation 0, LJ Oi OH llj were interpreted using 
a theory in which the tunneling dynamics was controlled 
by dynamic hyperfine and dipolar fields |l2|. A key con- 
nection between theory and experiment is the existence 
of 'hole-digging' in the distribution M(£,t) of magneti- 
sation M over longitudinal bias energy £ in the low T 
regime where the molecules behave as 2-level systems 
and phonons are unimportant 

Th e hole 

appears because a molecule can then only tunnel when 
in resonance. Without nuclear spins resonance only oc- 
curs if the local energy bias |£| < A Q . However the nu- 
clear spins mediate inelastic tunneling over a much larger 
range |f | < £ D ; in Fe-8 at low fields, £ /A ~ 10 5 . The- 
ory thus predicted an initial holewidth £ , provided £ <C 
the range of energy bias caused by the dipolar field dis- 
tribution. The experimental half- width w of the hole at 
short times, in initially strongly annealed samples, was 
interpreted 0, ^| to be the total half-width E Q of the 
multiplet of nuclear spin states coupled to the molecule- 
this latter can be calculated if the hyperfine couplings are 
known 

8]- 

Support for this interpretation came when iso- 
topes were varied ; w agreed both with approximate 
calculations of E Q , for different sets of nuclear isotopes, 
and with values of £ Q extracted from independent mea- 
surements of the sample relaxation rate Tq 1 . 

It is important to understand this hole-digging, not 
only to evaluate the claim that it demonstrates the role 
of nuclear spins in controlling the molecular dynamics 



(and the use of the holes as diagnostic tools), but also 
because the physics of the hole formation turns out to be 
rather interesting. There are 4 obvious questions which 
one may ask: 

(i) Why is £ = E Q — w in the _Fe-8 experiments, 
when M(0) is small? 

(ii) The experimental hole width w(t) varies in time 
(with w(t) — > w at short times) - how and why? 

(iii) w(t) also depends very strongly on the initial po- 
larisation M(t = 0) in the system (M(t) = f d£M(£,t)), 
but w is apparently independent of M(0) for small 
M(0)/M s ; is this expected theoretically? 

(iv) Why is the experimental lineshape Lorentzian? 

In fact most published experiments are on partially an- 
nealed samples, where M(0)/M s < 1 {M s being the sat- 
urated magnetisation of a fully polarised sample- hence- 
forth we measure M(t) in units of M s ). In the theory 
@ [l1 the quantity £ D was not necessarily the same 
as E ; and neither the hole lineshape, nor the depen- 
dence of w(t) on M(0) or t, were given. Moreover the 
theory only deals with weakly or strongly polarised sam- 
ples, whereas almost all experiments are on partially an- 
nealed samples, with intermediate polarisation- these are 
beyond the reach of analytic work. 

To address these 4 questions, we present here numerical 
results for the hole lineshape, and the holewidth w(t) as 
a function of t and initial annealing. We also calculate 
results for the f e-8 system, to provide predictions which 
can be compared with experiment. 

1. Short-time hole width: The experimental interpre- 
tation, that at short times w = E Q , appears to be based 
on the assumption (made explicitly in ref. |lj|) that 
the nuclear bias potential behaves as a noise, fluctuating 



2 



rapidly over the whole nuclear multiplet (halfwidth E Q ). 
Actually this assumption is unnecessary- in many cases 
(lif the hole width depends on the energy exchanged be- 
tween nuclear spins and the molecules during inelastic 
tunneling If the strength of the hyperfine inter- 

actions between the molecular spin S and the fc-th nu- 



clear spin Ik 



k ' 



then El = EM) 2 (h + 1)4/3. 



However energy is exchanged over a range £ Q ~ nui , 

where oj = N ~ 1 Ylk u; k * s the mean hyperfine coupling 
to the N nuclear spins in the molecule. In low fields 

2k = XX^fc" /^k) 2 j where = 9k I^nHo is the Zeeman 
coupling to a field H Q . In experiments with zero applied 
field, H a is just the strength of the local dipolar field 
from the other molecules- its root mean square in nearly 
annealed Fe-8 molecules is ~ 0.05 T. If k > v~N, then 
M to -> E . 
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FIG. 1: The variation of k(H ) with transverse field H^, 
calculated for a single Fe-8 molecule at low fields, for various 
isotopic mixtures. The hyperfine interactions were assumed 
to be dipolar, except for the core polarisation coupling to 
57 Fe, which was taken from experiment. 

Evaluation of n at low fields is simple for Fe-8 if we 
assume that all hyperfine couplings are dipolar except 
those to 57 Fe nuclei ^lj, since all nuclear positions and 
moments are known IT?! . The couplings {wj!}, {oJ^} can 
then be computed 18j. Fig. 1 shows the results, for 
different isotopic concentrations in Fe-8. We see that if 
it were not for the strong intermolecular dipolar fields, 
would be much less than F Q ; however, in fields ~ 0.05 T, 
they will be very close to each other. 

2. Time dependence of holewidth: To find the time 
dependence of w(t) we perform Monte Carlo simulations 
for the time evolution of M(£,i) on a simple cubic lat- 
tice, averaging over several hundred random initial spin 
configurations consistent with the initial M( t — 0) (for 
previous MC simulations see refs. 0, 0] (polarised 
systems) and (strongly annealed systems)). Indi- 
vidual spins, modelled as point dipoles oriented along 



±2, relax incoherently in a local bias field £ at a rate 
r^ 1 ^) = To^e-I^ , where r" 1 cx Al/A N , with A the 
tunneling matrix element and A^r a nuclear spin diffu- 
sion rate 0|- The dipolar fields at each site, calculated 
numerically, are updated at time intervals St -C t (typ- 
ically St/ t ~ 10~ 3 at the beginning of the decay) by 
flipping individual spins with probability 1 — e~ St / TN ^\ 
We define in the usual way the magnetisation distribu- 
tion M(£,t) = J d 3 r[Pi (£, r; t) - F x (£, r; t)] for the sam- 
ple, where P CT (£,r;t) is the probability that a molecule 
at position r with spin orientation a =| , J, will be in bias 
field £ at time t. A typical result is shown in Fig. 2. 
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FIG. 2: The tunneling hole dynamics, obtained by Monte 
Carlo simulations, performed on a model cubic lattice of 40 3 
point dipoles {En is the strength of the nearest neighbour 
dipolar interaction, and Ed/^o here). The main graph shows 
the time evolution of M(£, t), defined in the text, for a sam- 
ple in external field H z = 30£ o at M(0) = 0.25 (for which 
M(£,0) is nearly Gaussian, with a halfwidth Wd ~ 4.4Fd 
for this cubic lattice). The two upper insets are (i) at left, 
the magnetization curve M{i) vs the yt and (ii) at right, the 
tunneling hole half-width w(t) as a function of time. The in- 
set at the bottom right enlarges the region around £ = to 
show the hole dynamics. Note that time t is shown in units 
of t and bias £ in units of £ . 

A number of important general features of these results 
are as follows. If and only if the Gaussian halfwidth Wjj 
of M(£, t = 0) satisfies Wd ^> £o, we find: 

(i) for t < tq the hole halfwidth w(t) — > w , and w a — 
£o identically. Given the result above, that £ D ~ E a , this 
means the short-time holewidth w n is indeed measuring 
E Q , as previously assumed jirilllll20l|. 

(ii) When t > tq the hole rapidly broadens. The sharp 
rise in w(t) when t ~ t d appears not to be an artifact of 
the simulation. We have not yet been able to establish 
the asymptotic form of w(t) for t/t ^> 1, but for 1 < 
i/r < 5, w(t) grows slower than t 0A . 

(iii) For t > t q a square root relaxation (ie., \M(t) — 
M(0)| cx {t/TQ) 1 / 2 ) is always found, for any M(0). This 
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is obeyed until the fractional change in \M(t) — M(0)| ~ 
0.1-0.15. 

3. Hole lineshape: When M(0) < 0.5, we find the 
hole lineshape is Lorentzian. If M(0) > 0.5, deviations 
begin to appear; when M(0) ~ 1 the lineshape is not 
Lorentzian at all. 

It is important to understand the origin of this 
Lorentzian lineshape. In a strongly annealed sample the 
dipolar field fluctuates very rapidly from one molecule to 
another, but on a coarse-grained scale the sample should 
look homogeneous. We can therefore define the average 
time Tde(S.) for a molecule with initial bias £ to come to 
resonance. The total magnetization is then given by the 
homogeneous ensemble average: 

M(t) = M eq + (M(0) - M eq ) J deiV(0e-*/ T "(« (1) 

where M eq = M(t -> oo) and N(£) = J d 3 r[P T (£, r; 0) + 
Pj_(£,r;0)] is the total dipolar field distribution. The 
Lorentzian hole lineshape results directly from a 
Lorentzian form for tZ (£). We demonstrate this by cal- 
culating Tde(0 using the Monte Carlo procedure. 
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FIG. 3: The average decay rate tJ r (£) obtained by the MC 
simulations in a cubic sample of 40 3 molecules (diamonds). 
By solid line we plot the fitting Lorentzian curve with the half 
width £o- Three insets are: (i) the magnetization relaxation 
curve M(t) vs the vt; (ii) the average decay time T de (£) vs 
the £ 2 (diamonds); and (iii) the set of points in space, where 
for M(0) = 0.15 at t = molecules are in resonance. 

At t = all molecules were split into different bias 
groups of width £n/2. When some molecule flipped, the 
number of molecules in its group was reduced by one. 
The time taken for each bias group to decrease its initial 
number of molecules by a factor e is assumed to be the 
decay time r^ e for this group. Fig. 3 shows the results 
for M(0) = 0.15; one has a Lorentzian line-shape with 
half- width £n. 

Note one simple implication of this result: when t/ro > 
1, the dominant contribution to the integral Q comes 



from £ > £o where Td e (Q ~ £ 2 ; by making the integral in 
Q dimensionless, we see the magnetisation must decay 
like \M(t) - M (0)| ~ Thus the Lorentzian lineshape 
of the hole in M(£), in strongly annealed samples, is di- 
rectly connected to the square root relaxation, arising 
from the dipolar form of the intermolecular interactions. 
Notice, incidentally, that the distribution of molecules 
in the sample found to be in resonance for these strong 
annealings is indeed almost homogeneous (Fig. 3, inset 
bottom left). 

4- Application to Fe-8 system: Instead of point dipoles 
we now have 8 spin-5/2 Fe +3 ions correctly positioned 
and oriented inside each molecular unit cell, arranged in 
a triclinic lattice array. The lattice parameters, easy and 
hard axis orientations, were taken from standard sources 
|l6l,ll7j . As an example we show in Fig. 4 results for M(t) 
and t ( ^, 1 (^), for an initial magnetisation M(0) = 0.2. 
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FIG. 4: Plot of t^ 1 , as in Fig. 3, but now for a crystal 
of Fe-8 molecules (with 40 3 molecules arranged in a tri- 
clinic lattice). The inset above shows the decay of M(t) 
with \fi for a partially annealed sample (initial polarisation 
M(0) = 0.2). The lattice parameters for the Fe-8 triclinic 
lattice (a=10.522(7)A, b=14.05(l)A, c=15.(l)A, with angles 
a = 89.90(6),/? = 109.65(5), 7 = 109.27(6), in degrees), and 
the positions of the 8 spin-5/2 ions in each unit cell, were 
taken from ref. |l7|. The bias energy £ is now given in Kelvin 
units, and the value £ D = 6 mK is that appropriate to natural 
isotopic concentration [Tlj . 

Our basic conclusion here is that none of the results 
found for simple cubic lattices with point dipoles is essen- 
tially altered by going to a magnetisation distributed over 
the 8 Fe ions, or by the change of lattice structure- one 
still has a Lorentzian lineshape for r^, 1 , and the \ftpr~Q 
relaxation form for M (t), for M(0) < 0.5, for the simu- 
lations we have done [22j. 

5. Nearly polarised samples: Finally we examine the 
opposite limit of strong initial polarisation, which has al- 
ready been discussed both analytically and numerically 
|12ll2l| . In a sample of generic shape, the molecules that 
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can tunnel inelastically (ie., those satisfying the 'reso- 
nance' condition |£| < £ D ) form well-defined 'resonance 
surfaces' in space if the sample is nearly polarised, and 
£ <C Ed |2lLl23| . This is because now the dipolar field is 
dominated not by short-range fluctuations in space, but 
by the smoothly varying demagnetisation field. These 
resonance surfaces can be visualised by plotting those 
molecules satisfying the resonance condition at time t, 
during a Monte Carlo simulation. The spatial variation 
of the longitudinal component of the dipolar field can 
then be described by a slowly-varying demagnetisation 
field E(r,t), plus a smaller rapidly va rying: fl uctuation 
5t;(r,t), having Lorentzian distribution |12L l2- r | over val- 
ues of 5£, and {8£(r,t)) = 0. 
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FIG. 5: The main figure shows the time evolution of the dis- 
tribution M(£, t) of magnetisation caused by dipolar fields, in 
an initially polarized sample (M(t — 0) = 1); the same cubic 
sample of Fig. 2 was used. The 3 insets show: (i) the reso- 
nance surface at t = 0, showing all molecules in the sample 
satisfying the resonance condition in the text; (ii) the same 
result at t = 4r , where we see that the molecules satisfying 
the resonance condition have now spread in space (iii) the re- 
laxation of the normalized magnetization M(t) vs VI for the 
same sample (with the square root form shown as a thin line) . 

It then becomes interesting to ask how these surfaces 
evolve in time. Again we must resort to Monte Carlo sim- 
ulations. Fig. 5 shows a typical example- for all sample 
geometries we have studied the surfaces break up rapidly 
once the fraction (I — M(t)) ~ 0.07 — 0.11, correspond- 
ing closely with the onset of deviations from the square 
root scaling. This is intriguing because it indicates the 
deviations may be caused not only by higher correlations 
between the molecules, as previously assumed 0; this 
result needs further investigation. Notice that the hole 
lineshape is now irrelevant to the square root form of 
the relaxation dynamics, which comes instead from the 
Lorentzian distribution of <$£(r, t) 0,E^|. The hole line- 
shape is now controlled by the sample shape. 

Our essential conclusion in this paper is that one can 



use Monte Carlo methods to investigate the relaxation 
characteristics in regimes inaccessible to analytic theory. 
This allows one to understand the lineshape, width and 
time evolution of the hole that appears in M(£, t). 
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